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Abstract 

Fractional systems with Riemann-Liouville derivatives are considered. The 
initial memory value problem is posed and studied. We obtain explicit steer- 
ing laws with respect to the values of the fractional integrals of the state 
variables. The Gramian is generalized and steering functions between mem- 
ory values are characterized. 
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^] ', 1. Introduction 

O I Several approaches to the generalization of the notion of differentiation 

to non-integer orders are available in the literature of fractional calculus. 
Those include, e.g., Riemann-Liouville, Griinwald-Letnikow, Caputo, and 
the generalized functions approach. Mathematical fundamentals of fractional 



calculus are given in the monographs [la . |2Q| . l2ll l25l l26l . |29| . The idea of 



fractional calculus is used by engineers to model various processes since late 



1960's 11, 27, 30, 32, 33 



Here we study fractional-order systems using a vector state space repre- 
sentation, where we include also initialization. Our systems are described by 
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initialized Riemann-Liouville derivatives Jsj9| . Fractional state space repre- 
sentations of systems are given, e.g., in |l|, y|, |9|, |l2|, |3l|. Particularly in [9|, 
the method of initialized fractional calculus is addressed |6l M. IqL Il5 . 

The initial condition problem for fractional linear systems is a subject 



under strong consideration [14|, [15|, |16|, [17|, l22|, |23|, |24|, |28|. Several authors 



claim that the Riemann-Liouville derivative leads to initial conditions with- 
out physical meaning. Others, e.g. Heymans and Podlubny [10], contradict 
such claim, giving several examples of physical meaning to the Riemann- 
Liouville initial conditions, by introducing the concept of "inseparable twins" . 
Proper initialization is crucial to the solution and understanding of fractional 
differential equations. An important tool in the solution of fractional differ- 
ential or differintegral equations is provided by the Laplace transform. We 
consider the problem of solving initialized fractional differential systems us- 
ing the Laplace transform. This can be interpreted as the initial memory of 
the system. Then we develop the idea by looking at a given system not in 
the perspective of the state but also throughout memory, that can be of a 
different fractional order than the order of the system. 

The paper is organized as follows. Section [2] gives a brief review of frac- 
tional integrals, the fractional Riemann-Liouville derivative of scalar func- 
tions, properties for special functions in fractional calculus, as the Mittag- 
Leffler function and its generalizations, the matrix case, and the a-exponential 
matrix. Then, in Section [31 we formulate the fractional system E^ under 
consideration, where a denotes the fractional order of the system and ip 
the initialization function. We prove the form of solution for such kind of 
systems, we introduce the notion of /3-controllability, and we give a gen- 
eralization and description of the controllability Gramian. Our notion of 
/3-controllability is strongly connected with the system and the given ini- 
tialization function ip. The new possibilities are, essentially, present in the 
special case when 1 — a = (3. Indeed, for 1 — a = (3 = we obtain the classical 
linear dynamics but we are comparing the integral of the state. Moreover, 
for integer-order systems (a = 1) the initialization term in the Laplace trans- 
form is the Laplace transform of a constant and becomes ^(s) = a/s; for 
fractional-order systems this term is time- varying into the future, that is, the 
history of x(t) has the interference into the infinite future [9j. The classical 



notion of controllability [34| is obtained from our definition by taking a = 1 
and (3 = 0, what is interpreted as controllability without memory. Our no- 
tion is, however, more than a generalization to fractional-order systems: it 
introduces a new meaning, where we do not observe the state of the system 



but the integral of state along time. Finally, in Section 2] we state the no- 
tion of steering control and prove its optimality. We end with Section \5\ of 
conclusions and some future perspectives. 

2. Preliminaries 

In this section we make a review of notations, definitions, and some pre- 
liminary facts which are useful for the paper. Since several different notations 
and various definitions are used for fractional derivatives, our goal here is to 
make precise the definitions and formulas for solutions of linear differential 
fractional systems with Riemann-Liouville derivatives and initial conditions 
given by the value of a fractional integral at to- We recall definitions of frac- 
tional integrals of arbitrary order, the Riemann-Liouville derivative of order 
a G (0, 1), and a description of special functions in the fractional calculus. 
We present also the idea of initialization response with connection to the 
initialized linear fractional differential equations for order a G (0, 1). As in 
the classical (non-fractional) linear case, the Laplace transform provides a 
strong tool. We give here some useful formulas connected with the matter. 



Definition 2.1 ( [13|, |26|, [29| ) . Let ip G L x ([t , h],R). The integrals 
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I? Mt) = —l virKt-rr-'dr, t > t , 



r a) 



to 
h 



I£_^) = — -/ <p(T)(T-t) a -'dT, t<t 



't 

where T is the gamma function and a > 0, are called, respectively, the left- 
sided and the right-sided fractional integrals of order a. Additionally, we 
define l9 + = I® 1 _ := I (identity operator). 

We have the following formula of integration by parts for fractional inte- 
grals. 

Proposition 2.1 ([13]). Let a > and 1/p + \/q < 1 + a, p > 1, q > 1, 

with p 7^ 1 and q ^ 1 in the case 1/p + 1/q = 1 + a. For (p G L p ([t , ti], R) 
and ip G L q ([t , t±], M.) it holds: 



¥>(T)Jg + V>(r)dT = / ^{r)I? x Mr)dr. (1) 

to J to 



Remark 2.1 (|2j). When t Q = we write Iq + = I a and then I a f(t) = 
(f*f a )(t) ; where ip a (t) = f^y fort > 0, <p a (t) = fort < 0, and tp a -)> 5(t) 
as a — > 0, with S the delta Dirac pseudo function. 



Definition 2.2 (|13|, |26[). Let ip be defined on the interval [to,ti]. The left- 
sided Riemann-Liouville derivative of order a and the lower limit to are de- 
fined through the following: 

D i^ = r^{i)' 'j[V)(«-')— '*• 

where n is a natural number satisfying n = [a] + 1 with [a] denoting the 
integer part of a. Similarly, the right-sided Riemann-Liouville derivative of 
order a and the upper limit t\ are defined by 

Remark 2.2. Let a G (0,1). For functions <p> given in the interval [to,tx], 
each of the following expressions 

D °^ = r<lhj| { «** ~ ^" iT = It ( W» ■ 

°--^ = -fo^jt r ^ - ^ dT = -i t c£w» 

are called a fractional derivative of order a, left-sided and right-sided respec- 
tively. 



From 13|, Theorem 2.4] we have the following properties: 

Proposition 2.2. If a > 0, then D^ o+ I" Q+ ip(t) = <p(t) for any ip G Li(t ,ti), 
while It + Df + ip(t) = (p(t) is satisfied for ip G I£ )+ (L 1 (t ,t 1 )) with 

J t Q 0+ (L 1 (t ,t 1 )) = Mt): ?(*) = ig + V(t), V»eL!(tb,ti)}. 
However, 

In particular, for a G (0, 1] we have 



Similar results hold for right-sided fractional derivatives. There exist also 
a formula for the composition of Riemann-Liouville derivatives (cf. [26| ) . 



The next proposition is based on 29|, Corollary 2, p. 46] and is particularly 
useful for our purposes. 

Proposition 2.3. The formula 

[ * f{t)D? Q+ g(t)dt = f ' g(t)D°_f{t)dt , < a < 1 , (2) 

J to J t 

is valid under the assumption that f(t) G i| * (L p ) and g(t) e I? + (L q ) with 
l/p+l/q<l + a. 

Proof. The equality (J2J) follows from ([I]) if we denote Df f(t) = f(t), 
Dt o+ g{t) = ip(t), and take into account that I^ + D^ + f(i) = fit) is valid 
for f(t) E I« + {L{) (cf. [29]). ' ' D 

Another important result is the formula for the Laplace transformation of 
the derivative of a function ip (cf. |4J). Like the Laplace transformation of an 
integer order derivative, it is easy to show that the Laplace transformation 
of a fractional order Riemann-Liouville derivative is given by 



ra-l 

£ [D^(t)] {s) = s a C [<p{t)\ (s) - J2 

b— n 



ra-l 

k 



CkS k , n — 1 < a < n , 

fc=0 

where c k = (/o+^VWW For a e i * 1 ]* we have C [^o+V^W] (s) = 
s a £ [y(i)] (s) — a, where according with many authors the condition a = 
(/o-f °V(£)) t _ is considered a technical initial condition, mainly due to the 
Laplace transform and without a good physical interpretation p, |7|, |8|, |9_|, [l4j, 



15|, [161, 0, 122|, |23|, m \2M 



In papers of Hartley and Lorenzo p, |7|, |8|, |9[ , one can find the fundamentals 
about the idea of initialization process. By d^cp we denote the derivative 
without initialization or just the Liouville derivative defined by the inverse 
Laplace transform: 

dMt) ■= £- 1 [s a F(s)] (t) , 

where F(s) = C [<p(t)] (s). For a fixed function (p and a > 0, by initialization 
function ip = ip(t) we mean a function that satisfies the following: 

%(^„^)(<)+^(t). 



Here we consider only problems starting at 0. Then we have 

For a function x : [0, T] — > M. n , we use similar notation as in the classical 



case: 



D% + x(t) = D, 



o+ 



/ Xl (t) \ ( D% +Xl {t) \ 

x 2 {t) D« + x 2 (t) 

\ Dg + x n (t) J 



\ x n (t) J 

Such situation, when for each component we use the same fractional order 
a of differentiation, is called in the literature the fractional-order derivative 
with commensurate order. Commensurate order is used in connection with 
both Riemann-Liouville or the Caputo derivative (cf. |3l. Il2l. |31|). 

We need to mention now two important functions and their extensions to 
matrices. 

Definition 2.3 ([5]). The two-parameter Mittag-Leffler type function is de- 
fined by the series expansion 

a > 0, 0> 0. 



Let A e 



E a ,p(z) = J2 
By 



fe=0 



T(ka + /3) 



E a ^(At a ) = J2^ 



4-ka 



fc=0 



T(ka + /3) 



it is denoted the extension of the two-parameter Mittag-Leffler function to 
matrices. 

For — 1 we obtain the Mittag-Leffler function of one parameter: 

„jfc 



EJz) 



£ 



z 



T(ka + 1) 

In the matrical case, E a (At a ) = YlT=o ^ k r<ka+i) ■ 
Definition 2.4 ([13]). Let A e R nxn . Then, 

t ka 



ei % = f 



-^ 



fc=0 



T[(k + l)a] 



» t (k + i)*-i 

E A 

fc=0 



T[(k + l)a] 



t a - L E a , a (At a ) 



denote the a-exponential matrix function. 



For a = 1 both functions given in Definitions 12.31 and 12.41 are equal and 
coincide with the classical exponential matrix: E a (At a ) = e^* = exp(Ai). In 



19j the following properties are proved: 



Proposition 2.4. Let a > 0. Then, 

(a) Dfre? = Ae£; 

(b) D%_<p(t) = Acp(t), where <p(t) = S(T - t) = ei {T ' t] ; 

(c) I + f*Ae^dT = E Q (At a ); 

(d) f t E Q (At a ) = Aei\ t > 0; 

(e) D? o+ (l* + x(t))=x(t). 

Additionally, we need a more general form for the item b) of Proposi- 
tion [231 In fact this item is also true for S(t) = if**?- 1 E^ a+fi {At a ) . 

It is easy to explain the following properties of fractional integrals of the 
a-exponential function. If to = we write I a = i? =0+ . 



Q-l 



(a) Let k G NU {0}. Then, I^^Si - r(ka+1) , 

(b) p-"e? = E a {At); 

(c) I^e^ = t a ^- l E a>a+fi (At a ); 

(d) I?e M = I^eiU = t^E ltfi+1 (At); 

( e ) (/ /3 e^) 4=0 = 0,/3>l-a; 

(f) (i^ef a) t=Q = a,ae R n . 

Let S(t) = t a+ P~ 1 E 0itCt+ p(At a ). The next result follows from the technical 



lemmas proved in [19 



Lemma 2.5. Let < a < 1 and p,(t) e Iq{L^), where l/p+ 1/q <l + a for 
p such that all components of S(T — t) belong to I^_{L p ). Then, 

f S(T-r)D^fi(r)dr= f AS(T - r)fi(r)dr . 
Jo Jo 



3. Systems and the Controllability Gramian 

We consider an initialized fractional linear time-invariant control system 

Dg + x(t) = Ax(t) + Bu(t) , £>0, (S Q> ^) 



') n X rn 



where a E (0, 1], x{t) E W 1 , u{t) E M m , matrices A E R nxn and B E 
together with an initialization vector ip(t), t < 0. We consider piecewise 



constant controls «(■)• In |15| authors explain that the "state" of the system is 
not given by the dynamic variable vector. This is because of the initialization 
vector. The derivative notation is 

D% + x(t)= d? X (t)+m, 

where dfx(t) is the uninitialized fractional derivative starting at t — and 
ip is the initialization vector function, determined as 

ip(t) = lim — — / ; v \ dr 

Comparing with |9>], we have a = — oo: we are involving the history from the 
left axis. In the particular case when x(t) = a for — oo < t < 0, we obtain 
the following: 

, . . ,. d ( 1 f° a , \ at 

m = lim - — / -dr 



oo dt \r(l-a)J u (t-r) a J r(l-a) 

Proposition 3.1. Let < a < 1. The forward trajectory of the system E Qi ^ 
evaluated at t > has the following form: 

T M,«)= / S(t-T)(Bu(T)-i>(r))dT, 
Jo 

where S(t) = e^*. Moreover, for x{t) = a, t E (— oo,0], we have: 

^{t,a,u) = E a {At a )a+ / S(t - T)Bu{r)dr . (3) 

Jo 



Proof. First part is proved in [9| using the Laplace transform. We only need 
to prove here the second part. We have 

X(s) = (Is a - A)- l BU(s) - (Is a - Ay^is) . 



Applying the inverse Laplace transform for 

-(Is a - A) _1 *(s) = (Is a - A)- 1 s a - 1 a 

we get 

C' 1 [-(Is a - A)- 1 *(s)] (t) = E a (At a )a . 

Hence the formula ([3]) holds. □ 

We introduce now the more general situation when one has some external 
measure of the memory. This is represented by the fractional integral of order 

(3. 

Definition 3.1. Let (3 > 1 — a. The memory of order (3 of the forward 
trajectory j(t,i(),u) evaluated at time t > is defined by 

M fi {t,i/>,u) := ig+iit^u) . 

Additionally, we define the value of the memory of order (3 at t = as 

Proposition 3.2. Let (3 > 1 — a. Then for the system E aj ^ we have 

Mpit, ij),u)= f ®p{t - r) {Bu{r) - ^(r)) dr for t > 0, 
Jo 



where 

Mt) = t a+p - l E ata+ ^At a ) = f+P- 1 J2 A 



k l 



Moreover, for ip(t) — a for t G (— oo, 0] we have 

Mp(t,ip = a,u) = t? E a< p(At a )a + / $g(£ - r)Bu(r)dr . 

Jo 

Proof. Using the Laplace transformation we get 



X(s 



oo 



C [M (t, i>, u)] (s) = ^ = - Y, A k s^ k+1 ^-^(s)+J2 s^ k+1 ^ p A k BU(s)) 



s>- 

k=0 k=0 



The desired formula follows from the fact that 



Mt) = £ _1 



k -(k+l)a-p 



5> 

,k=0 

oo 

t a + P-lJ2A h 
k=0 

t a+ P- l E a , a+p {At a ). 



(t) 



ika 



T((k + l)a + /3) 



Moreover, for \&(s) 



s a ~ 1 a = C 



r(l-a) 



[s) we get that 



£ 



-i 



J- A fe s -(fc + l)a-/3^ 



k=0 



/r 1 



(*) 



5> 



fc „-(fc+l)a-/3 



a 



_k=0 

t^E a ^At a ). 



(t) = t"J2 Ak 



ika 



k=0 



T(ka + /3 



a 



Definition 3.2. Let T > 0. T7ie system S a ^ zs controllable with memory 
of order (3 > 1 — a on [0,T] if there exists a control u defined on [0, T] such 
that 

Mp{t^,u)\ t=T = b. 

We denote by 

Qt= [ (T - ^-"-^(T -t)BB*$*p(T -t)dt 
Jo 

the /3- controllability Gramian on the time interval [0,T] corresponding to 
the system £ Qj ^. As in the classical case, Qt is symmetric and nonnegative 
definite. 

Theorem 3.3. Let T > and Qt be nonsingular. Then, 

(a) for b e M. n the control law 

u(t) = -(T - t)^- a -VB*^(T - t)g T 1 /r(^ ! b) , (4) 



where fx{ip,b) = —b — J Q 3>p(T — T)ifj(r)dT, drives the system to b 



in 



time T. 
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(b) Among all possible controls from L 2 a ([0,T],R m ) driving the system to 
b in time T, the control u defined by 0) minimizes the modified energy 
integral 

£^{u) := / \(T - t) a+ P- l u(t)\ 2 dt . 
Jo 

Moreover, 

S^(u)= f \(T-tr^- l u(t)\ 2 dt=<Q T 1 f T (^,b)J T (^b)>, 
Jo 

where < ■, ■ > denotes the inner product, and 

f T (ij,b) = -b- [ $ (T-T)i>(T)dT. 
JO 

Proof. From the form of u we directly have that 
Mp(T,i;,u) = f T (^b) + b 

-U(T- t) 2(1 - Q ^^(T - t)BB*^(T - t)dt\ Q^Mi?, b) 

= f T {il>,b) + b-Q T QT 1 Mi/>,b) 

= b. 

Item (a) is proved. Similarly to the classical situation [34], to prove item (b) 
we begin noticing that 

i-T 

/ \(r - ty+^u^ 2 dt 

Jo 

= f I (T - t) l - a ^B*^{T - t)Q^f T ^, b)\ 2 dt 
Jo 

= f \ T _ t \2ix- a -P) (B*^(t-t)Q^f T (iP,b),B*^(T -t)QT X h^b))dt 
Jo 

= If \T -t\^- a -^(T -t)BB*%{T -t)dt,Q^f T (TP,b)\ 

= (QTQ T 1 fT^,b),Q T 1 f T (^b)) 
= (f T (^,b),Q^f T ^,b)) . 



11 



T 


T 



Let us take another control u for which (T — t) a+l3 1 u(t) is square integrable 
on [0, T] and Mp(T,if),u) = b. Then, 

- / (T - i) 2 ^- 1 ) < M (t), (T - ^-"-^b^t - t)g T V T (^, 6)) d* 

u(t),B*$i(T -t)Q?f T (iJ>,b))dt 

o 

Hence, 

(T_t)2(i-«-« („(<), u( t ))dt = / (T-t) 2{1 - a -^(u(t),u(t))dt 
o </o 

and we obtain 

-T 

(T-t) 2 ^- a ~^\u{t)\ 2 dt 

o 

JO 

which gives the intended minimality property for the integral. □ 

Example 3.1. Let £ a =o.5,y> be the following system in M?: 

' D° f Xl (t) = x 2 (t) , 
D° fx 2 (t) = u(t) . 

Let us take 6=1 n ) ■ Since ^4=1 n n ) > -^ = ( ) > an d ^ — 0, we 

obtain for t = the formula for the solution with the initialization ip in the 
following form: 

l(t,if>,u)= i | («(r) - ip(r)) dr . 

JO \ y/n(t-T) 



12 



Let [5 >0.5. Then, 

t f)-0.6 t P 



*e(t) 



r(o.5+/3) rfi+/3) 







f 



-0.! 



T(0.5+/3) 

For simplicity let us take /3 = 0.5. Then, 

$0=0.5(0 = ( "f 

and t/ie Gramian is the following: 

( 2T 2 AT 1 - 5 
4T 1 - 5 j. 

3v^F 

TTius, n(£) = p /T^t , 1 Qj* 1 fT^ify. Let us take the end point b equal to the 



r(3/2) • 



zero vector and ^(0 — 0. T/ien ; n(£) = I gX^r-t ) Jo V'i( r )^ T ; an <^ 

£/ie norm 



T \ 2 /.T 



withu{t) = ^ (Jo V'iCO^' 7 ") ■ 

4. Steering Laws 

In the classical theory it is quite easy to explain why the rank condition 
of matrix B equal to n is sufficient for controllability. In the case of fractional 
order systems we need to use a special function. 

Proposition 4.1. Let rankB = n, matrix B + be such that BB + = I, 
and gp(-) be the matrix function such that $p(t)gp(t) = I, t > 0, and 
lim <&p(t)gp(t) = I . Then the control 



u(t) = T{a + ^ 



B+g(T ~t)(b + J ^(T - r)^r)dr\ , t 6 [0, T] , 



T 
transfers system S a ^ to b in time T > 

13 



Proof. It follows by direct calculation: 

>y(T,il>,u) = - [ ^/s(T-r)^(r)dr 
Jo 

+ ^j <$> p (T-t)BB + g p (T-t)L + j ^^T-T)^r)dr\dt = b. 

D 
Example 4.1. Let us consider the system 

E^ : D% + x(t) = u(t) , 
where D* + x{t) = d% + x{t) + ip(t). LetbeR,T>0. Then, 

ftx+p-1 

*^ = WWY B+ = B = ^ 

i <- T 



M(T, ^ u) = / (T - t)**- 1 (u(t) - m) dt . 

i[a + p) Jo 



The steering law 

u(t) = r( ^ /3) (T - ty-^f (b + £ <MT - r)i,{r)dr 

transfers the system to b. 

If the rank condition is satisfied, the control u given by (@| drives the 
system to b at time T. Our goal now is to find another formula for the steering 
control by using the matrix [A\B] instead of the controllability matrix Qt- 
It is a classical result that if rank [A\B] = n then there exists a matrix 
K G M(mn,n) such that [A\B] K = I G M(n,n) or, equivalently, there are 
matrices K\, K 2 , . . . K n G M(m, n) such that 

BK X + ABK 2 + ■■■ + A n ~ 1 BK n = I . 



In [19j we introduced a convenient notation for composition of fractional 
derivatives with the same order a, where a G (0, 1]. Let us recall here the 
notation for fractional Riemann-Liouville derivatives. Let i?Q^_/x(t) = //(£). 
We define, recursively, R^ +1 /j(t) := D« + (R^^(t)),j G N. 



14 



Theorem 4.1. Let (3 > and rank [A\B] = rank [B,AB, . . . , A n ~ x B\ = n 
for system E^. Let p be such that QpiT — t) G I^_(L p ) and let tp be a real 
function given on [0, T] such that 

(1) £<p{t)dt=l; 

(2) RotKt) e Jo+(A?) M j = 0, . . . , n - 1, where 

tit) = 9p(t) (b + J *p{T- r)ij(r)dr] <p{t) , 

and $^(T - t)gp(t) = I n for t G [0, T], while l/p +l/q<l + a. 
Then, the control 

u(t) = K lt i(t) + K 2 D« +f i(t) + ■■■ + KnESp^tit), t G [0, T] , 

transfers the system E ai </, to the value b at time T > 0. 

Proof. Using j — 1 times the formula d2J) of integration by parts and Lemma 12751 
we get for j = 1, . . . , n that 

/ <$> p (T-t)BK j FQ>f 1 ii{t)dt= f $0(T -t)A j - l BK^{t)dt. 
Jo Jo 

Then, 

/ $ p (T-t)Bu(t)dt= / & p {T-t)Bfj,{t)dt 
Jo Jo 

and finally 

7 (r,V,«) = - / ^(T-r)V(r)dr 



.7 
■T 



+ / $ m {t)[b + j $f,{T-T)il>{T)dr\<p(t)dt = b. 



a 



5. Conclusions 



In this work we develop the theoretical aspects regarding optimality and 
steering laws for fractional systems E Q ^ with Riemman-Liouville derivatives 
of order a and a given initialization function ip. The novelty of the theory 

15 



here introduced is the "memory state point of view" for studying fractional 
systems. Main results provide explicit steering laws that drive the system 
from the initial memory value characterized by ip to a desired memory state 
of the system. Moreover, we obtained the optimal law in the sense of mini- 
mization of the modified energy. 

The present work opens several possibilities of future research. Interest- 
ing open questions include the study of memory observability conditions for 
initialized systems and realization issues. The realization question is par- 
ticularly interesting. Roughly speaking, after solving a fractional system of 
order a we neglect that we know the true state, hiding it under an additional 
integral of a different fractional order. From the point of view of transfer 
functions, it means that we do not know of what type, in the sense of the 
fractional order and class of derivative, a realization can be done. One pos- 
sible approach is to realize the transfer function using not only the direct 
dynamic but coming back to the dynamic of the additional integral of the 
state. This research is in its beginning, and it will be developed in the future. 
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